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Abstract 

Effect of hydrophilicity of the confining walls on lamellar phases in oil- water- 
surfactant mixtures is studied in a slit geometry. In contrast to strongly 
hydrophilic or hydrophobic walls, which induce parallel orientation of lamel- 
lae, the lamellae can be oriented perpendicularly to the neutral walls when the 
material properties and the ter mo dynamical state of the sample are suitably 
chosen. When the elastic energy associated with compression or decompres- 
sion of the lamellae parallel to very weakly hydrophilic walls is sufficiently 
large, then changes of the film thickness lead to a switch from the parallel 
to the perpendicular orientation of the lamellae. Our general arguments are 
confirmed by explicit mean-field calculations in a lattice vector model. 
In self-assembling systems such as diblock copolymer melts, binary or ternary surfactant 
mixtures, lipids etc., various ordered structures occur spontaneously on a nanometer length 
scale Jl|-§[. The most common is the lamellar phase [|J. In ternary surfactant mixtures the 
lamellar phase consists of oil- and water-rich layers separated by monolayers of surfactant. 
It is natural to expect that the layers are parallel to the confining walls and the parallel 
orientation of lamellae is commonly assumed. Indeed, near a strongly hydrophilic or hy- 
drophobic wall a water-rich or an oil-rich layer is adsorbed at the surface respectively, by 
which the orientation of the subsequent layers is fixed |5|-f|. However, when the surface is 
neutral, i.e. it is neither hydrophilic nor hydrophobic, then neither component is favored 
and no orientation of amphiphiles is preferred near the wall. Therefore perpendicular rather 
than parallel orientation of the lamellar phase may be stable between two parallel surfaces 
that are neither hydrophilic nor hydrophobic (see Fig.l). In the case of very weakly hy- 
drophilic walls perpendicular orientation was indeed observed in Monte Carlo simulations 
H of the phenomenological Landau- Ginzburg model ||. Perpendicular orientation should 
be favoured when the distance between the walls and the period of the lamellar phase do 
not match. On the other hand, when the distance between the lamellae fits the distance 
between the weakly hydrophilic walls, parallel orientation can be stable. 

Films with the lamellae oriented parallelly or perpendicularly to the walls have different 
mechanical, electrical, optical etc. properties. Therefore the possibility of inducing a switch 
between the two orientations by changing a control parameter such as the distance between 
the walls L is potentially very important for various applications. For example, when some 
electrolyte is added, the electrical conductivity between the confining walls is high in the 
perpendicular orientation, whereas the parallel orientation resembles a sequence of capaci- 
tors. The switch of the lamellae orientation leads also to an abrupt change of the refraction 
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index and the reflectivity, and of elastic properties of the film. A lamellar film containing a 
few lamellae has the thickness L ~ 30 4- lOOnm. Such a lamellar film would correspond to 
a switch of a mesoscopic size - a device which may be useful for modern nanotechnologies. 
In this letter we show under what circumstances such a switch from the parallel to the 
perpendicular orientation of the lamellae can occur. 

The confined fluid induces the solvation force between the walls |10| . 



where Q,p,A,T and \i are the grand thermodynamic potential, the bulk pressure, the area 
of the confining walls, the temperature and the chemical potential respectively. / can 
be directly measured in surface force apparatus experiments 14]. In simple fluids 



oscillating / reflects packing effects of molecules for L up to several molecular diameters, 
since the streched or shrunk confined fluid tends to assume the bulk structure. Similarly, 
but on a much larger length scale, / reflects packing effects of oil- and water-rich domains in 
the complex fluids. For L > 3A, where A is the period in the bulk phase, lamellae are parallel 
to the strongly hydrophilic (or hydrophobic) walls f||7]||. When L and A do not match, the 
lamellar structure is stretched for L > or shrunk for L < Ln where Ln = NX + A/2 
denotes the equilibrium thickness, corresponding to no stress, for a given number N of 
periods of the confined phase. For L ^ the confined lamellar phase responds elastically 
to the stress strain, just like joined strings, each string representing one period of the lamellar 
phase j7|. The elastic contribution to the excess thermodynamic potential 

Vt ex A = O - uj b LA, (2) 

where u>b is the bulk thermodynamic-potential density, has a form 

n% = BAL 2 /2L N , (3) 

with AL = L — L N . The elastic modulus B depends on the amphiphilicity of the surfactant 
and on the thermodynamic state ||. 

Because neither component is favoured near the neutral wall, it has a disordering effect 
on the lamellar order due to the missing neighbors JT5J . The missing-neighbors contribution 



to the excess potential Vt ex depends on the orientation of the lamellar phase (see Fig.l). For 
the parallel orientation (||), Q ex = the interactions are cut along a single layer or along 
an interface. For the perpendicular orientation (_L), Q ex = Q ex , the water- water, surfactant- 
surfactant and oil-oil interactions are missing at the surface in proportions depending on 
A and the thickness of the monolayer a. The density distributions near the surfaces tend 
to minimize the effect of the missing neighbors at the walls, and determine the values of 
the surface tensions mi and <jj_. Sufficiently far from phase transitions the effect of the wall 
should extend to distances comparable with the characteristic length of the system, which 
in this case corresponds to A. If the lamellae are perpendicular to the confining walls, there 
is no constrain on the period of the lamellar phase, which can be equal to the bulk value, 
and in this case there is no elastic contribution to Q e f. Thus except for narrow slits, i.e. for 
L > 3A, we expect Q e f ss 2a± independently of L. In the case of the parallel orientation of 
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lamellae, for L ^ Ln, there are deformations of the structure in the whole slit like in the 
case of the hydrophilic walls and m 2a\\ + Q^f (L). 

By suitable tuning of hydrophilicity of surfaces, amphiphilicity, T and surfactant volume- 
fraction p s one could in principle obtain comparable cry and a±. For weakly hydrophilic walls 
one can expect <j\\ < cr±. Thus, the parallel orientation should be stable in the case of no 
stress (L = Ln). The perpendicular orientation may be stable for slightly different wall 
separations, corresponding to stretched or shrunk parallel lamellar phases, if the elastic 
contribution to QjP is sufficiently large. According to the phenomenological theory [TJ 



B = 97T 2 (kT) 2 A/2[64k(A/2 — a) ] , where k is the bending elasticity of monolayers. Thus 
the compression or decompression induced switch to the perpendicular orientation is to be 
expected for relatively small A, rather than for highly swollen lamellar phases. 

Whether the switch can indeed occur in particular systems should be verified either 
experimentally or by a more accurate analysis. The theoretical analysis should allow for a 
determination of cr±,<y\\ and Q e i (or more accurately Q?, x , Q e f) in particular systems. There 
are several theoretical approaches to bulk self-asembling systems. The approach suitable 
for a particular problem depends on the relevant length scale. In our case two length ratios 
are important, namely L/X and A/a. L/X is the number of oil (water) layers in the slit 
for the || orientation, and A/a characterizes the liquid-crystalline lamellar phases. For A/a 
sufficiently large ( A/a > 4), the water- rich layers are sufficiently thick to be compressible 
as in bulk liquid. Experimental results [P,[7|,pj}-|I^| show that the lamellar phases respond 
to compression or decompression by shrinking or swelling the water-rich (and/or oil-rich) 
layers, rather than by developing deformations characteristic for solids or by deformations 
of the molecular structure of amphiphiles. Hence the degrees of freedom related to the 
structure of molecules can be neglected, as long as L/a > 2. Dislocations, frustration etc. 
can be expected for stiff, solid-like structures, i.e. when the water-rich (oil-rich) layers are 
very thin, A/a < 4. Since for our problem the relevant length rations are 1 < L/X < 10 and 
4 < A/a < 10, the molecular degrees of freedom, frustrations etc. can be disregarded. On the 
other hand, the popular approaches in which the monolayers are approximated by infinitely 
thin mathematical surfaces may be oversimplified, particularly for comparable thicknesses 
of the layers of water and surfactant (4 < A/a < 6). 

In a presence of a wall the surface term describing the interaction of the surface with 
the water, oil, and amphiphiles in various orientations must be included. In a purely phe- 
nomenological approach one cannot apriori know what values of additional surface param- 
eters are physical. Therefore one cannot know if the predicted phenomena, occurying for 
certain values of parameters, are physical or not. It is thus advantegeous to consider a 
semi-microscopic approach and relate the interactions with the surface to the effective in- 
teractions in the bulk. For example, for the water-covered surface the interactions between 
various components and the surface are determined by the interactions between these com- 
ponents and water and there is much less ambiguity in defining the surface contribution in 
such an approach. 

All the above discussed requirements for the appropriate description of confined am- 
phiphilic systems are fulfilled by the semi-microscopic CHS (Ciach, H0ye and Stell) lattice 
model ||17|| . The model is described in detail in Ref. |I7| and here we remind it only briefly. 
The microscopic states are /3;(x) = 1(0) if the site x is (is not) occupied by the state i, where 
i = 1, 2, 2 + M denotes water, oil and surfactant in different orientations respectively. In 
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the case of close-packing and oil - water symmetry only one chemical potential variable is 
relevant, namely \i = — fi sur f = li% — fi S urf, with // snr / = /ij for i > 2. The Hamiltonian 
in a presence of external fields /ij(x) can be written as: 

tf = ^EEM x )u^(x-x')p,(x') + 

+ £ £ hi (x)p,(x) - ft £ (ft (x) + p 2 (x)) . (4) 

x j x 

The lattice constant a = 1 is identified with the length of the amphiphiles (~ 2nm). Nearest- 
neighbor interactions are assumed and U#(x — x') vanishes for |x — x'| ^ 1. In the case 
of oil- water symmetry the water- water and oil-oil interactions are of the same strength — b, 
and the water-oil interaction energy is set to zero. The interaction between an amphiphile 
in an orientation w at x and a water (or oil) particle at x' is —cuo ■ Ax (or +cu) • Ax), where 
Ax = x — x', thus opposite orientations of amphiphiles are preferred by the water and the 
oil particles, as in real systems. Two amphiphiles with orientations u> and &' at x and x' 
respectively contribute —g(u x Ax) • (a)' x Ax) to the system energy (when |Ax| = 1), 
i.e. g supports formation of planar monolayers with amphiphiles parallel to each other and 
perpendicular to the surface they occupy. Explicit expressions for Uy (x — x') defined above 
can be found in Ref . §J . For lamellar phases the direction perpendicular to the planar layers 
of surfactant, oil and water is distinguished. We arbitrarily choose the orientation of the 
unit vector normal to the layers and denote it by n. Then all orientations of amphiphiles 
such that uj ■ n > (or u - n < 0) are mapped onto a single state — ► (or <— ). Hence, once 
the direction of density oscillations is fixed, there are 4 states corresponding to water, oil, 
— > and <— at every lattice site. 

The model parameters characterize the effective interactions between different compo- 
nents, and the interactions with the surface can be expressed in terms of these parameters. 
When the surface is water-covered, all the components interact with the wall in the same 
way as with the bulk water, and /ij(x) can be easily derived from the bulk interactions. 
When the hydrophilicity of the wall is decreased, the interactions with the wall can be de- 
creased by the same factor h s for all the components. Thus, < h s < 1 should correspond, 
at least qualitatively, to physical surfaces, ranging from neutral through weakly hydrophilic 
to strongly hydrophilic walls. 

The model parameters are not directly measurable, but the calculated values for mea- 
surable quantities can be compared with experiments. For example, the model parameters 
can be choosen so that A/a is the same as in a particular experimental system. In this way 
the model parameters can be related to different experimental systems in the bulk. It turns 
out that once the bulk properties in the model and experiments are the same, for strongly 
hydrophilic walls (the case with no free surface parameters) the modulus of compressibility 
B (in kT units) in the model and experiments agree very well ||. Therefore, in contrast 
to purely phenomenological models with arbitrary parameters, the CHS model can quite 
reliably predict phenomena not yet studied experimentally. 

The structure of the stable confined lamellar phase in mean-field approximation (MF) 
corresponds to the global minimum of the grand-thermodynamic potential 

n MF (T, fi, l) = E M x ) ( fcT ln (p*( x )) + + fti(x) - tiki + • (5) 
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Here 0«(x) = Z) x / Z)j Uy (x — x')pj(x') is the mean field and p«(x) is the MF-average of 
Pi(x). In the slit geometry we assume 1 < z = £3 < L. We consider two classes of density 
distributions, corresponding to the parallel and the perpendicular orientations of the lamellar 
phase. In the first case the two distinguished orientations of amphiphiles are Cj = ±2. In 
the second case we assume that the layers are parallel to the (y, z)-plane, and Cj = ±x. The 
slit is infinite in directions x = x\ and y = x-i- In the case of the parallel orientation of the 
lamellar phase p«(x) = Pi(z). In the second case, due to the deformations of the structure 
near the surfaces, the density distributions are Pi(x) = Pi(x, z). In the direction x we assume 
periodic boundary conditions in a system of a size A, since Q ex for such system is the same as 
for infinite periodic structure with the period A (we also considered smaller and larger sizes 
to be sure we find the stable structure). We use the method of finding local minima of Q MF 
by solving numerically a set of self-consistent equations for densities, tested in earlier works 
0. Then we find the structures corresponding to the lowest values of fijj* and Q e f, and 
eventually we calculate — Q ex to find the global minimum. It turns out that Q^ x — Vt ex 
depends sensitively on the hydrophilicity of the surfaces h s , on the ratio of water-surfactant 
c/b and surfactant-surfactant g/b interactions and on the period of the lamellar phase, which 
for given interactions depends on p (or p s ) and T. We find that for highly swollen lamellar 
phases the parallel orientation is preferred, and that for smaller periods A the perpendicular 
orientation stabilizes between neutral walls. We find that there exist systems for which the 
switch of the orientation does take place. Fig. 2 refers to such a system and corresponds to 
the stability region of the lamellar phase far from phase boundaries, with A/a = 6. The 
switch takes place for relatively small A and for weakly hydrophilic walls, as expected. Note 
that in the case of the neutral walls the period of fijj^L) is A/2 (Fig. 2a), since either oil- or 
water-rich layers can be adsorbed at the surfaces. In contrast, in the case of the hydrophilic 
walls only the water- rich layers are formed near the walls and the period of Q^ X (L) is A (see 
Fig.2c). 

In MF the correlations between the lamellae are neglected. Such correlations are par- 
ticularly important close to the melting of the lamellar phases where they lead to creation 
of channels connecting the neighboring water-rich layers and eventually to the transition to 
the microemulsion. Far from phase boundaries, however, the correlations should play less 
important role. In the case of the hydrophilic walls the MF results for the CHS model agree 
quite well with the results of experiments. The effect of the correlations should not depend 
significantly on the kind of the confining walls and we expect that the correlations can have 
quantitative, but not qualitative effect on the results. 

Our results show that the switch is to be expected for A/a ~ 6, i.e. with the water (oil) 
layers twice as thick as the surfactant monolayers, and sufficiently far from phase boundaries. 
From the rough relation 2a/A ~ p s we obtain the estimation p s « 0.33±0.1 for the surfactant 
volume fraction which should correspond to the occurence of the switch. For this range of 
p s the lamellar phase is stable at room temperatures for example for extensively studied 
experimental system such as water, decane and C W E 5 |L8| . 



Various experimental methods can be used for detecting the switch, due to the abrupt 



change of various propertis of the sample. For example, SFA [TT[ measurements should 
show vanishing force for the perpendicular orientation and elastic responce for the parallel 
orientation. Reflectivity, index of refraction and electrical conductivity change abruptly 
when the switch takes place. Measurememts of the above quantities should give significantly 
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different results for the two orientations, thus allowing for observation of the phenomenon. 
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Fig.l. Schematic illustration of parallel (a) and perpendicular (b) orientations of the 
lamellar phase in a slit. Shaded regions represent oil-rich domains, while white regions 
represent water-rich layers and thick black lines represent the surfactant monolayers. 
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Fig. 2. Q^ X (L) (white circles) and Q e f(L) (triangles down) as a function of the wall 
separation (in units of a) for c/b = 2.4, kT/b = 2.8, /i/b = 3 and g/b = 0.15. (a) h s = 
(neutral walls); (b) h s = 0.015 (weakly hydrophilic walls); (c) h s — 1 (strongly hydrophilic 
walls) . 
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